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^-H . Abstract 

o 

^^ ' We consider the version of the Zee model where both Higgs doublets 

^D ' couple to leptons. Within this framework we study charged Higgs 



Q^' decays. We focus on a model with minimal number of parameters 

consistent with experimental neutrino data. Using constraints from 
neutrino physics we (i) discuss the reconstruction of the parameter 






^ . space of the model using the leptonic decay patterns of both of the 

5t ! two charged Higgses, h^2 ~^ ^t^i^ and the decay of the heavier charged 

Higgs, /i^ — > hfh^; (ii) show that the decay rate T{hf — > ^x^Vi) in 
general is enhanced in comparision to the standard two Higgs doublet 
model while in some regions of parameter space T{K^ -^ fJ-^i^i) even 
dominates over T{hf — > r'^i'i). 



1 Introduction 

Neutrino oscillation experiments, including the results of KamLAND P] have 
confirmed the LMA-MSW oscillation solution of the solar neutrino problem. 
Together with the earlier discoveries in atmospheric neutrinos 2J, one can be 
fairly confident that all neutrino fiavours mix and that at least two non-zero 
neutrino masses exist. 

In the standard model neutrinos are massless. Among all the existing 
models to generate small neutrino Majorana masses the seesaw mechanism 
j3] is perhaps the most popular. However, this is not the only theoretical 
approach to neutrino masses. Other possibilities include Higgs triplets jl], 
supersymmetric models with broken i?-parity [SlEl; some hybrid mechanisms 
that combine the triplet and the i?-parity ideas ^ and radiative mechanisms 

iHiin!. 

Here we consider a particular radiative mechanism, the Zee model . In 
this model the scalar sector of the standard model is enlarged to include a 
charged SU(2) gauge singlet scalar and a second Higgs doublet. This particle 
content allows to write an explicit lepton number (L) violating term in the 
scalar potential and leads to neutrino masses at one loop order. In the 
Minimal Zee Model (MZM), only one Higgs doublet couples to leptons [Til] . 
As a result, dangerous Flavour Changing Neutral Current (FCNC) processes 
are forbidden. It has been shown ^I] that combining SNO, KamLAND and 
K2K experimental data this version is ruled out. 

However, this does not mean that the Zee model is ruled out. The original 
version, from now on called the General Zee Model (GZM) ^2], in which 
both of the two Higgs doublets couple to the matter fields has been shown 
fn\ IT^ IT^ to be consistent with atmospheric and solar neutrino data as well 



Once one allows both of the Higgs doublets to couple to leptons the 
number of model parameters increases. Here instead of working with all the 
couplings of the model we will consider a scheme, previously discussed in 
references fn\ I13j. where the neutrino mass matrix has a two-zero-texture. 
This particular GZM will be called Next to MZM (NMZM). 

In the Higgs sector, after spontaneous breaking of the electroweak symme- 
try, the charged gauge singlet mixes with the charged components of the two 
Higgs doublets. The resulting charged Higgs eigenstates {hf with i = 1,2) 
decay to states with charged leptons and neutrinos. These decays can be 
used, in principle, to reconstruct the Majorana neutrino mass matrix. 

We will show that due to the constraints imposed by neutrino physics, 
the Br{h^ -^ Yli ^if^~^) is enhanced in comparision to the two-Higgs doublet 



models (2HDM) of type-I and type-II ^. Moreover, we will show that in large 
parts of the parameter space Br{hf —>■ ^^ i^ifJ^'^) ^ Br{hf -^ J2i ^i'T^)- For 
details see section IHl 

The rest of this paper is organized as follows. In section |21 we give the 
generalities of the GZM and work out the Higgs mass spectrum of the model. 
In section El we study charged Higgs production at a future e"'"e~ collider. In 
section m we discuss the bounds on the parameters of the model coming from 
FCNC processes constraints. In section El we describe the Majorana neutrino 
mass matrix within the GZM and in the NMZM. In section El we discuss the 
connection between neutrino physics and charged Higgs decays. In sectional 
we present our conclusions and summarize our results. 

2 The Model 

2.1 Generalities 

If no new fermions are added to the standard model neutrino masses must 
be always of Majorana type, i.e. the mass term must violate L. In the Zee 
model an L = 2 charged scalar, /i"*", is introduced. Since this field carries 
electric charge its vacumm expectation value (vev) must vanish. Therefore 
in this model L cannot be spontaneously broken. However, h^ can be used 
to drive the lepton number breaking from the leptonic sector to the scalar 
sector. In order to accomplish this a new SU(2)2, doublet has to be added, 
as a result an explicit L violation term can be written. This term is given 
by 

liec^pH^H^h- + H.c. (1) 

where yU is a coupling with dimension of mass and Hi and H2 are doublets 
with hypercharge Yi = ¥2 = 1. 

The most general Yukawa couplings of the model can be written as 

-Cy = Hll,),,H,eRj + e„^L°/,, C{L'')^h- + H.c. , (2) 

where Lj are lepton doublets, crj are lepton singlets, C is the charge conju- 
gation operator, H^ (a = 1, 2) and / are 3x3 matrices in flavour space, e^^ 
(a,/3 = 1,2) is the SU(2)l antisymmetric tensor and i,j = 1,2,3 are family 
indices. / is an antisymmetric matrix due to Fermi statistics. 

^In type-I only one of the Higgs fields couples to the SM fermions, in type-II one Higgs 
field couples to up-type quarks and the other Higgs field couples to down-type quarks. 
There is another version called type-Ill jSl where both Higgs fields couple to all SM 
fermions. 
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In general both of the two Higgs doublets can acquire vev's, (Ha) = v, 
with V = ^Jv\ + v\ ~ 246 GeV. As usual, the ratio of these vev's can be 
parametrized as tan/? = v^jvx. 

1.1 Higgs potential and scalar mass spectrum 

Though in this work we are interested mainly in the charged Higgs sector of 
the model and its relation with neutrino physics, we will briefly discuss the 
full scalar mass spectrum. 

The Higgs potential is invariant under a global SO (2) transformation 

E'\ ^ / cos/5 sin/3\ (EA 

Moreover, the Yukawa Lagrangian given in Eq. (0) is also invariant under 
the above transformation if the Yukawa matrices are appropriately rotated, 
namely 

'H'A _ / cos/? sin/3\ /Hi' 



U'J \^-sin/3 cospJ\U2j' ^^^ 

Therefore the full model is invariant under global SO (2) transformations of 
the two Higgs doublets. Thus, we are free to redefine our two doublet scalar 
fields by making an arbitrary SO (2) transformation. A particular choice of 
fields corresponds to a choice of basis. There is a basis in which only one of 
the two Higgs doublets acquire a vev. In the context of the 2HDM of type-Ill 
it is called the Higgs basis. Notice that in this basis tan/? = 0. 

In the Higgs basis the most general gauge invariant scalar potential of the 
model, consistent with renormalizability reads 

V = ^ilH\H^+^xlHlH2-[^llH\H2 + YL.c.] + ]^\^{HlH^f 

+h2{HlH2f + \^{hIh^){hIh2) + \{hIh2){hIh,) 

+ l]^\,{H\H2f + [\{H\H^) + \j{hIH2)]H\H2 + H.c. 

+Aio|/i+|'(/f|/f2 + H.c.) + fieo^pH'^H^h-. (5) 

Since we will not deal with CP-violating effects we only consider real coeffi- 
cients 



Minimization of the scalar potential, Eq. (jSj), leads to the conditions P^ 

2 -'- \ 2 

^4 = 2^6^^^- (^) 

These conditions can be used to ehminate /i^ and ^\ as independent variables 
from V . 

Of the original ten scalar degrees of freedom, three Goldstone bosons 
{G^ and G°) are absorbed by the W^ and Z^. The remaining seven physical 
Higgs particles are: two CP-even (/i° and H^ with m/io < ttiho), one CP-odd 
(A°) and two charged Higgs pairs (hf and hf). 

In the basis $^ = (G", -ff~, h~) the squared-mass matrix for the charged 
Higgs states is given by 



^^ = I M^± -f,v/V2 I , (7) 



where 



^0 -fiv/V2 MI3 



Ml± = /i2 + 2^'A3 



-M33 = ^A + v'X,. (8) 

The matrix element M|^± corresponds to the squared-mass of the charged 
scalars {H^) that in the absence of the S\]{2)l singlets k^ would be physical 
Higgs particles. 

The squared-mass matrix M.'jj can be diagonalized by the rotation matrix 

/I \ 

i? = I COS V9 sin (y9 . (9) 

yO — sin (/? cos (pj 

where the angle y? characterize the size of the H^ — h^ mixing. 

The mass eigenstate basis in the charged Higgs sector is defined as H"'' = 
(G^, h~[ , /i^) and the rotation angle is given by 

\/2vjj 

Here Mi and M2 stand for the masses of the scalars h^ and /ig which are 
given by 

Ml, = 1 (^M^, + Ml, T ^(M^, - Ml,y + 2/^2^2^ . (11) 



«in2(^ = 77^^- (10) 
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Figure 1: Production cross section for charged scalars h^ at an 1 TeV e~^e~ 
collider with unpolarized beams. 

In the Higgs basis the squared- masses for the CP-odd and CP-even Higgs 
states are given by [12] 



M 



AO 



^HO,hO 



1 



:^'(A5 - A4) 



Mlo + v\X, + As) ± J[M% + v^{X, - Ai)]2 + Av^Xl 



(12) 



3 Charged Scalar Phenomenology 

3.1 Cross section 

In the following we discuss charged scalar hf {k = 1, 2) production at a future 
e^e~ collider, /i^ are produced in e^e~ annihilation via s-channel exchange 
of a 7 or Z^ ^. The total cross section for the process e'^e" -^ h^h^ will be 
the sum of three terms 



cr,. 



(Tj + (TZ + (TjZ 



(13) 



^There is also a t-channel Yukawa production through neutrino exchange but due to 
the smahness of this contribution to the total production cross section we do not consider 
it here. 



corresponding to the pure photon, pure Z and photon-Z interference contri- 
butions respectively. Thus 

(15) 



is - M2)2 + Mlrl 

1 

1927r 



z"- z 

4„2 



3 9 ^w nur^ o„2^^ i , a. .2 
s-Ml, 



^z, = -T7^(i''-r{Wt,-2si){-l + As 



'-Z 



{s-MlY + MlTl 
where s^, = sin6'^, c^ = cos6'u,, 



(16) 



/lf2 

/3=yi-4-^ (17) 

and Wii = W22 = cos ip and Wu = —W21 = sin ip. 

From Eqs. (jH)), (fTT|) and (fT2j) it can be noted that fixing Mi 2 does not 
fix M^o/jO and M^o. Therefore, it is possible to take Mi 2 and Wik as free 
parameters without being in confiict with LEP bounds for the CP-even and 
CP-odd Higgs masses |IZl CHI CHI • 

In Fig. n we show the cross section at an 1 TeV e~^e~ collider with unpo- 
larized beams. There we have taken 80 Gev < M^ < 500 GeV. The spread in 
the plot is due to the dependence of the cross sections on ip. It is important 
to notice that for small (large) values of (p the cross section for h^ increases 
(decreases) while the cross section for /i^ decreases (increases). Figure Q il- 
lustrates the situation for the case ip = 0. In that case h^ coincides with 
the SU(2) doublet H+ (solid line) and ht with the SU(2) singlet h+ (dashed 
line). The dotted-dashed line corresponds to cosip ~ 0.6. 

In Fig.[T]it can be seen that up to a mass of ~ 350 GeV the charged scalars 
have a cross section larger than 10 fb. Assuming an integrated luminosity of 
1 ab"^ this implies that at least 10*^ charged scalar pairs will be produced. 

3.2 Decay Widths 

After the spontaneous electroweak symmetry breaking charged leptons ac- 
quire mass, namely 

M, = -^^t;„n,. (18) 

* a 
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In the mass eigenstate basis for the charged scalars we have 

— C-Y 3 z/ij'Oj/je/jj (cos (/?/;, j['~ — sin y^/ig") 

-^{yuf C{2fij)eLj{^\w^hX + cos^/i^) + H.c. (19) 

where, in general, the couphngs Oij are given by 

^ ^tan/? — 1 „ , , 

O = -V2 ^M, + -n2. 20 

V cosp 

Charged scalars /i^2 "^ill decay through the couplings Oij and fij. Possible 
leptonic final states are i^dt- Possible final states involving quarks are diUj. 
These decays are determined by the couplings O^^ where, in general 

O^ = -Vt^M, + ^lil (21) 

V cosp 

Here q refers to up-type and down-type quarks, Mq are the diagonal quark 
mass matrices, and 112 ^^^ 3x3 Yukawa coupling matrices of the second 
Higgs doublet. Notice that in the Higgs basis M^ = (f/-\/2)ni and O = I\.2- 
We are interested in the widths and branching ratios for leptonic fi- 
nal states. The Lagrangian (fT^ determines the two body decays hi2 -^ 
(Sj ^i)^^ ■ The decay rate reads 

nK - (E ^^n) = ^ Et^J^'^ + (2^)^w^i.] • (22) 

i i 

The couplings h^W' Z and h'^W~'~f do not exist in the Zee model. This can 
be understood as follows: since h'^ is a mixture of H~^ and /i"*" these couplings 
are determined by the SU(2) doublet component. However, in the 2HDM of 
type-III these vertices do not exist [20j. Therefore the decays h^ — > W~^j, 
W~^Z^ in the Zee model are not present at tree level. For this reason we do 
not consider them. 



4 Constraints from FCNC processes 

In the GZM FCNC interactions are induced by the charged and neutral 
Higgses. Bounds on the OjiOkm couplings can be obtained from the non- 
observation of tree-level processes i~ — ^ i'j'i'j^i:^. Constraints on OkiOkj 
come from radiative processes i~ -^ ijj induced by neutral Higgses. Limits 
on fikfkj and on Okifkj couphngs come from radiative processes mediated by 



Process 



Constraint 






e^e e 

' -^ e~fi~e~^ 

T~ — > e~7 



K 



|0i30n|< 1.3 X 10-3 (^)^ 
|Oi30i2|< 0.9 X 10-3 (^)' 
IO23O21I < 0.9 X 10-3 ^ *^" 



^ lOOGeV y 

1013021 + OsaOnI < 1.0 x lO'^ (^5^)' 

1013022 + O23O12I < 1.0 X 10-3 [j^^^ 

|Oi20n+02202l + 03203l| <4.1X 10-5/" 

|0i30n + O23O21 + O33O; 



xlO-^l" ^'° V 



•3303i|< 4.7 X 10-2 (^y 

/ » rf) \ 2 



IO13O12 + O23O22 + O33O32I < 3.3 X 10-2 



K 



X 



^ lOOGcV ) 



Table 1: Constraints on the parameters Oij from tree level and radiative 
FCNC processes induced by the neutral Higgs h^. 



charged scalars3. An important remark is that once the constraints on the 
fikfkj and OkiOkj couplings are satisfied the limits on Okifkj are no longer 
important, for this reason we do not list them. Tabled shows the constraints 
coming from the processes mediated by neutral scalars. Table E] summarize 
the limits on the fij parameters. Experimental constraints used in both 
tables were taken from 



5 Neutrino Physics 

5.1 Neutrino Mass Matrix in the GZM 

In this section we will discuss the neutrino mass matrix. The Majorana 
neutrino mass matrix in the Zee model arises at the one loop level through the 
exchange of the scalars h^ and /i2 as shown in Fig. |2l Assuming Mi, M2 ^ 

rrie, rrin, rrir we have 



iM,)ie = K[fij{Mt)jjOv, + Oij{M,)jjf, 



I n 



(23) 



■^These processes also give bounds on OikOkj- However they are weaker than those 
coming from radiative processes mediated by neutral Higgses. 



Process 


Constraint 


/i- -^ e-7 

r- -^ e-7 
r- ^ ;U-7 


/23/l3|<4.1xl0-^(,^)^ 

/23/i2|< 4.7 X 10-2 (,^)^ 
/i3/i2|< 3.3 X 10-2 (,,^)^ 



Table 2: Constraints on the parameters fij coming from radiative FCNC 
processes induced by the charged Higgs hf. 



where 



_ sin2v? fM^ 



(24) 



5.2 Neutrino Mass Matrix in the NMZM 

In this section we discuss the neutrino mass matrix in the context of the 
NMZM. In our scheme the neutrino mass matrix is assumed to be 



M„ 



K 



M,„ M 







e^ 



flT 



(25) 



M,^ M. 



flT 



The neutrino mass matrix, can be diagonalized by a matrix t/, which can 
be parametrized as 



U 



where Sij 




X 



Cl3 S13' 

1 

-Sl3 Cvi, 



C12 S12 0^ 
-S12 C12 
1, 



(26) 



sin6'jj and Cjj 



cos 9ij. Phases are zero since only real parame- 



ters are considered. 
From 



U^MM 



M,, 



(27) 



and taking the limit sin2 9i^ = 0, since experimental neutrino data require 
sin2 6'i3 to be small ^3] , we can find approximate analytical expressions for 
the atmospheric and solar mixing angle as well as for Aml^: 

2 



tan2 9 



23 



M„ 



M 



e/j 



tan 26* 



12 



^2"'^ 



M,. 



Mee + M, 



{2i 



IIT 



Aml^ 



K 



W 



(M,„ - M„). 



ht ,-'-. h 



' ^ 1 

Figure 2: Loop diagrams for Majorana neutrino mass. Here i = 1,2 

Due to our two-zero-texture mass matrix (Eq. ()25|) ) we have an inverted 
hierarchy neutrino mass spectrum |22] and therefore Mee — M^^. Thus the 
neutrino mass matrix, Eq. ()25p. has only three independent entries which, 
from Eqs. ()28|) . can be written in terms of tan^ 623, tan 2^12 and Aml^, namely 



Mee ^ M.r 



V^^h 



K tan 29 



12 



^J2Km 



Me>. ^ ,; , , jr , , (29) 



k(1 + tan^; 



23 j 



tan ^23 y2Am^ 
/t(l-htan^23) 

Assuming that there are no large hierarchies among the couplings On and 
Oij, terms proportional to rrie, in the neutrino mass matrix, (see Eq. ((221)) 
can be neglected. Thus we obtain Eq. (j^ with O23 = O32 = 0. Under this 
constraint the mass matrix depends on k and on the seven parameters 

/12, /l35 /23, O12, O13, O22, O33, (30) 



as can be seen from Eqs. (j23jl and ()25j) . By using equations in ()29j) we can 
write four of these parameters in terms of the other three. Equations (jH?)|) . 
(fHTj) . (|HHjl and (jH^jl . in the appendix, give the expressions for /12, /13, /23, 
Oi3, in terms of O12, O22, and O33. Note that both /23 and O33 must be 
different from zero. 

Next we will consider the cases for which Eqs. fj56|l . (jSZj), (J58|) and (J59|l 
can be expressed in terms of a single parameter. We will call these cases 
the one-parameter solutions. Since O33 cannot be zero (see Eq. (jSZj)) we will 
parametrize all our one-parameter solutions in terms of this coupling. This 
leaves us with only four possibilities: O12 = 0, O13 = 0, fu = 0, /13 = and 
the remaining parameters in Eq. (jHUj) different from zero in each case. We 
will show below that the first two lead to solutions with large O33, while the 
last two lead to solutions with small O33. 
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5.3 The one- parameter region 

The main point here is that in these two cases (small and large O33) not 
only neutrino physics but the decay patterns of hf are governed by a single 
parameter. This allows an analytical approach to the problem of identifying 
a particular collider signature that allows to distinguish between different 
regions in parameter space. In the following we will discuss the four possibil- 
ities mentioned previously and we will estimate the values of the parameters, 
consistent with neutrino physics as well as with FCNC constraints, in each 
case. This discussion will be useful in our analysis of the decays of hi pre- 
sented in section 16.21 

5.3.1 The large O33 case 

Choosing O12 = and O22 = ("^^/'^r)033, as in references jT21 UB]; Eqs. 
(J56|l . (jSZj), (J58|l and (J59|l are reduced to the one-parameter solution 



_[l + (2 + 4 tan22^i2) tan^aa+tan^^as] V'Ar^^r 1 6.3 x lO'^ 



2 v^ K tan2 26112 tan 6^23 (1 + tan 6'23) m^ O33 k O. 



33 



„ V2tan^23 x/Am^g 1 1.9x10-" ,^^, 

fi;(l + tan6'23) rUr O33 KO33 



f: 



23 



1 v^Am^g 1 1.2 X 10"" 



K tan 2^12 m^ O33 k O33 

l + tang23 ^ .„^ 

2 V2 tan 26^12 tan 6^23 

The last values in each equation are obtained using the best fit point value 
for each neutrino observable. 

An upper bound for k can be estimated using the fact that 

sin2^ (Ml\ ^ V2v^^ 1 ( Ml\ ^ V^W ^^^. 

(47r)2 \Ml ) (47r)2 M^ - Ml \mI ) ~ (47r)2 M^ ' ^ ' 

Therefore for M2 < 1000 GeV and |/i| < 500 GeV |22j, we have that \n\ < 
10-2. For example for Mi = 200 GeV, M2 = 300 GeV, and /z = 100 GeV, we 
have 

sin2y? = 0.7 and k = 3.6 x 10-^ (33) 

On the other hand, from the expression for /12 in Eq. (p?T|l and imposing 
/12 ^ IQ-^ a lower bound on k can be found. Choosing O33 < 10-^, we have 
that \k\ > 10-^. For example, for /i = 2 GeV and with Mi and M2 as in the 
previous case, we have 

sin2(^ = 0.014 and k = 7.2 x 10-^ (34) 
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Using the value of k given in Eq. ()33|) we have 

1.8 X lO-'^ ^ 5.2 X 10-9 ^ 3.2 x lO'^ 
/i2 ~ 7=; , /i3 ~ 7^ , 723 ~ 7^ • (35j 

t^SS t^SS t>'33 

Now instead of O12 = we choose O13 = 0. Again, as in the previous case, 
we take O22 = ('^/i/'"^r)033 and the best fit point values for each neutrino 
observable. With k given by (j^Hjl . Eqs. (J3B|) . (j37j) . (j^ and (j^ become 

1.5 X 10-6 5.2 X 10-9 3.2 x 10-^ 

J12 ~ 7^ , /l3 ~ 7^ , 723 ~ 7^ , t>'l3 ~ 0.02 Css, 

1^33 <-^33 ^33 

(36) 
which is basically the same result obtained in the case with O12 = (Eqs. ()35|)). 

From Eqs. (jHSj) and ()36|) it can be seen that all the parameters can be 
below 10-^, with a hierarchy of order 10^ between /12 and the others fij. In 
this way the constraints on the couplings coming from FCNC interactions 
(Tables HI and 121) are always satisfied. 

Note that /12 < 10-^ requires O33 > lO-''. Therefore the range of varia- 
tion of O33 is restricted to 10-^ < O33 < 10-^. For n small, as in Eq. ^^. 

O33 ~ 10-2. 

5.3.2 The small O33 case 

If we choose /13 = and, in order to define the one-parameter solution in 
this case^ O22 = 0, Eqs. (jSHl), (|SZ|), dSHI) and dMI) become 

(l + tan023) VAmig 1 3.6x10-^2 
/12 ^IT-^^ — \ TTT" n 7r~ ~ 7^ ' 



2 v^ K tan 29I2 tan ^23 ^r O33 k O 



33 



723 ~— ^7 7^ 7^ , 137) 

K tan 26'i2 m^ C33 k C33 

v^ tan 26*12 tan 6*23 m^^ „„^ 

t>'i2 ~ ,. , , ^—. U^i ~ 2/'G'33, 

(l+tan6'23) rra^ 

^ y2 tan 2^12 ^ _^ 

^13 ~~rT~;^ Q *^33 ~ l-0<^33- 

1 + tan 6^23 

The last values in each equation are obtained using the best fit point values 

for each neutrino observable. Note that O12 ^ 10-^ requires O33 ^ 4 x 10-^. 

On the other hand, from the expression for /23 in Eq. (jHTjl . if O33 < 

4 x 10"^ and we impose the bound /23 < 10-^ we have that k > 3 x 10-^. 

^This choice allow us to define the one-parameter solutions. However, we stress that 
our results does not depend on this choice. Our main conclusions hold for any O22 < O33. 
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Case 


O33 


K. 


/i (GeV) 


Large O33 
Small O33 


10"^ - 10-2 
10"^ - 4 X 10-3 


10-5 - 10-2 
3 X 10-6 - 10-2 


2-500 
0.2 - 500 



Table 3: Range of O33, k and /i for the one-parameter solutions in the NMZM 



For example, if we choose /x = 0.2 GeV, Mi = 200 GeV and M2 = 300 GeV 
we have 

sin 2v9= 1.4 X 10-3 and k = 7.2 x 10-^ (38 

For the value of k given in Eq. (jSHI), that satisfies the bound k > 3 x 10-^ 
with Ml = 200 GeV and M2 = 300 GeV we have 



/: 



1 X 10- 



12 



O 



f'. 



3.2 X 10-9 



23 



33 



O: 



(39) 



33 



Now instead of /13 = we choose /12 = 0. Using the best fit point values for 
each neutrino observable, O22 = and n given by Eq. (jH^ . Eqs. (j3Bj) . (jSTj), 
(j3Hj) and (j3^ become 



/: 



9.8 X 10 



-10 



13 



O: 



f: 



3.2 X 10-9 



23 



33 



o 



o 



12 



32O33, o 



'13 



1.6 



33- 



(40) 



33 



In both cases (/12 = or /13 = 0) we can have all the five parameters of 
order of 10-^ without any hierarchy among them. In fact, the case /12 = 0, 
considered here, is a particular case of the one studied in reference jTJ in 
which all the parameters Oij and fij are of the same order of magnitude. From 
Tables Q and |2l it can be seen that FCNC constraints are always satisfied. 

A lower bound on O33 can be obtained using the bound /23 < 10-^. 
Together with the upper bound estimated previously we have 10-^ < O33 < 
4 X 10-3. Notice that for smaller values of k, as the one in Eq. (jHH)) . the 
range of variation is more restricted, 10-^ < O33 ^ 4 x 10-^. 

It is worth noticing that there are no more possibilities in the one-parameter 
solution case. The large O33 case, obtained when either O12 or O13 are ne- 
glected implies a hierarchy among the non zero fij and, depending on the 
case, on O12 or O13. In the small O33 case, obtained when either /12 or /13 
are neglected, it is possible to have all the parameters at the level of 10-^. 
Table 01 shows the allowed range of variation for O33, k and /i in each case. 
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6 Neutrino and Collider Physics 

6.1 Determination of the Neutrino Mass Matrix Pa- 
rameters 

In this section we discuss how the charged scalar decays can give some hints 
about the parameters that determine neutrino masses and mixing angles. 
Charged Higgs decays are governed by the same parameters that control 
neutrino physics so, in principle, the information coming from these decays 
can be used to reconstruct the neutrino mass matrix. Outside of the one- 
parameter regions analysed in section El the number of parameters is large 
and since neutrino flavour cannot be determined the mass matrix cannot be, 
in general, reconstructed. Despite this, in the limiting case of small mixing 
(v? -C 1), the fact that the the mainly doublet state decays are dictated by 
the Oij and the mainly singlet state decays are controlled by the fij leads to 
a situation in which the reconstruction of part of the parameter space of the 
model is possible. 

The charged scalar singlet h~^ does not couples to quarks. Thus experi- 
mentally the mainly singlet state can be differentiated from the mainly dou- 
blet state by the fact that the branching ratio to final states with quarks 
{uidj) must be smaller for the former than for the latter. Our main assump- 
tion here is that all the decays that we are going to consider have a branching 
ratio in the order of at least per-mille. 

In the following discussion we will use the notation h'^^ for charged Hig- 
gses. Here d and s denote the mainly doublet and mainly singlet states 
respectively. Note that (i=l,s = 2or(i = 2,s = l are possible. Ratios 
of branching ratios for h^^ can be used to obtain information about the Oij 
and fij couplings. In the case of /it we have 



Br{hj^{Z,u,n) Y.^01 



2 



(41) 



(42) 



and for hf 

Brjht ^ (E. ^.W _ E. 4 
Br{ht -^ {Y.,i^,)it) EJL 

Corrections to both ratios are oc ip"^ <^ 1. The interesting point here is that 
despite the large number of parameters the relative size of the fij couplings 
can be obtained by suitable combinations of ratios of branching ratios, for 
example 

Br^ - Brl + Br^ ^ fj^ 

Br^-Brl + Brl /^ ^ ' 
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Figure 3: Ratio of branching ratios Br'^'^'^ / Br'^^'^ 



{Br 



Brl + 



Br^^)/{Br^ - Br^^ + Brl) versus A^a/Zla (left) and Br'^/Br^^ = Br{hj 



(E,^.)/i+)/i?r(/i+ -> {Ek^k)e^) versus 0^/0 
right). See text. 



(01, + 01,) /{Of, 



01 + 






For the Oij the situation is more 



with Sr/ denoting Br{hf — > (^j ^'j)^^ 

comphcated but even in this case some information can be obtained from the 

ratios of branching ratios. For example, the relation 



Br{hj-.{Z,'y^)^^^) 



O?, + 02 



22 



Br{hj 



[Ek ^k)e+) 01, + 01, + 01 



(44) 
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allows to determine the relative importance of the couplings involved in these 
decays. 

Figure El shows the ratios of branching ratios described above. Any devi- 
ation from the small mixing assumption would lead to a large dispersion. 

There are two limit cases of particular interest where the decays of h\^ 
are correlated with the neutrino mixing angles, On ^ O13 ^ O22 < O33 or 
0\z ^ 0\2 ^ O22 < O33. Figure m shows both cases. In the left plot the 
variables y, and x, are given by 



yi 



Xi 



Br^ - Brl + Brl 




Br'i + Brl 



m, 



m-r 



BriK ^ {J:,u,)t^ 



a/2 tan 26*12 



1 



tan 6*23 



(45) 
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Figure 4: Ratio of branching ratios indicated by the variables yi (left) and 
1/2 (right) versus the atmospheric and solar mixing angles indicated by the 
variables Xi (left) and X2 (right). See text 

In the right one the variables 1/2 and X2 are defined as 



2/2 



X2 



^BA 



Bri + Brl 



Br'i 
1 



Brl + Brl 






Br{hj ^ {Z,^k)fi+) 



72 tan 2^ 



'1 + tan% 



23 J 



(46) 
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Another important decay, if kinematically allowed, that could be used to 
obtain information about /i is /i^ — *> hfh^. The decay rate for this process 
reads 

r(/i+ ^ hth') = — ^ W 1 - 4^' 



Here 



and 



/^ 



1671 M2 

sin2y9 



Ml 



A = -^ sin a cos 2ip + v 

y/2 ^2 



(A22 - A33) 



(47) 
(48) 



A22 = A7 cos a — A3 sin a, 
A33 = Aio cos a — Xg sin a 



(49) 



where a is the mixing angle that define the two CP-even Higgs mass eigen- 
states, h^ and H^. 

Figure El shows the decay rate r(/i^ — > h^h^) versus /x^. There we have 
fixed M2 = 400 Gev, Mi = 150 Gev, M^o = 130 GeV and a = tt/G. /i is in 
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Figure 5: Decay rate Tiji^ -^ h^h^) versus /x^ for fixed values M2 = 400 Gev, 
Ml = 150 Gev and M^o = 130 GeV. 

tfie range 0.1 GeV — 8 Gev in order to ensure (p <^ 1. The dispersion is due to 
tfie presence of the other couphngs, present in the scalar potential (Eq. (J3J). 
Apart from allowing the approximate determination of /i, measurements of 
T{h2 — > h^h^) in the range indicated by Fig. El will indicate that the small 
mixing limit is realized. 

6.2 Hierarchy of charged Higgs leptonic decays 

In this section we will show that in the NMZM, the decay process h\ -^ 
(X]j ^i)f^~^ is enhanced in comparision to the 2HDM of type-I and type-II. 
Moreover, it is shown that in large parts of the parameter space, h^ -^ 
C^i ^i)f^~^ can be the dominant leptonic decay. 

In the one-parameter solutions, described in sec. El the parameters O12, 
Oi3, fi2, /i3 and /23 are governed by the parameter O33. In this way the 
Br{hf -^ {Y2i ^i)(^t) are functions of O33. In order to find expressions with 
no dependence on k or O33 and correlated with neutrino physics observables 
(tan 2^12, tan ^23) we consider ratios of branching ratios in the limits Oij 3> fij 



17 




(Z^i i^i)f^'^)/Br{hf -^ (^,. z/i)r+) as a function of O33 for 
2 obtained with 0.2 < /i < 500 GeV, Mi = 200 GeV and 



Figure 6: Br{h{ — 

3 X 10-*^ < K < 10 

M2 = 500 GeV. All the parameters fij and Oij satisfy the bounds shown in 

Tables HI and 121 For all the dark gray (green) points O12 < 10~^. See text 



and Oij ^ fij, namely 



Brjht ^ {Y.^ ^i)^^^) E J(Q»2 cos ^Y + (2/,2 sin ^f 
Br{ht 



(E. ^.)^+) E^ [(0.3 COS ¥.)2 + (2/,3 siu y,)^] 

'^ for Oij- > /y 



E.02 



J3 



2^. /.2 



(50) 



(51) 



for Oij < /j. 



Clearly from Eqs. (p?T|) or (jHTjl . Eq. (jKT|) depend only on the neutrino mixing 
angles and charged lepton masses. We will call the regions of parameter space 
with either Oij or /^ dominance correlation regions. Ratio of branching ratios 
in these regions are k independent (or /i independent). In general, outside the 
correlation regions, the independence on /i approximately holds, but there is 
a dependence on M2. 

Fig.inishows the ratio Br{hl — > (^^ Ui)n^)/ Br{hl — > (^^ ^i)T^) as func- 
tion of O33. For all curves, we have used the best fit point values for 
Am|3, and the solar and atmospheric mixing angles. We have taken also 
3 X 10-6 < K < 10-2 obtained when 0.2 < /i < 500 GeV, Mi = 200 GeV 
and M2 = 500 GeV. The correlation regions correspond to the flat parts of 
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the curves. The large O33 case determined by Eq. ^^ with Ou = cor- 
responds to the sohd hne in the right part of the plot, while the large O33 
case with O13 = correspond to the dashed line. In the same way, the 
small O33 case described by Eq. (p?7j) with /13 = corresponds to the dot- 
ted line in the left part of the plot, while the small O33 case with /12 = 
correspond to the dotted-dashed line. The scatter plot was obtained by 
searching for all solutions compatible with neutrino data at 3cr level, and 
keeping O22 = {m^/nir) O33. 

In the large O33 case described by Eq. (jSII), the correlation region for fij ^ 
Oij is excluded because the parameters fij are above the values consistent 
with FCNC constraints (see Tables ^ and |21)- For the other correlation region, 
for which 0^ ^ fij, we have 



^^^ I 

Br{ht ^ {Zi^^)r+) \m 



for Oij > fij. (52) 



As shown by the solid line at the right of Fig. IHl the contribution of fij can 
increase the ratio of branchings ratios up to a factor of 10^. In this way 
the decay hf —>■ (^^ i^i)fi'^ may become observable in future colliders. The 
dark gray (green) points were selected from the full scatter plot by choosing 
O12 < 10^^. They are well fitted by the solid line which represents the 
one-parameter solution with O12 = as given in Eq. ()31|) . 
In the small O33 case with /13 = 0, we have 

Br{ht ^ {E.u,)t-^) 

2 tan2 2^12 tan2 623 ml r ^ ^^ . 

Z 9 — T lO^^ Oij > Jij 

2tan2 2^i2 + (l + tan^23) < 



1 + 2 tan 023 + (1 + 8 tan^ 2^12) tan^ ^23 



8 tan^ 2^12 tan^ 9 



23 



for Oij < fi. 



f ITir 



2 



tan^ 023 ^ for Oij > /, 



\mJ ^'^ "-^ '-' (53) 



1 for 0„ « f 



ly 



As shown in the left part of Fig. IHl in this case the ratio of branching ratios is 
larger than one, and therefore an inverted hierarchy for the leptonic decays 
of the lightest charged Higgs is obtained. In this way, in the small O33 case, 
the most important leptonic decay channel for the charged Higgs h^ must 
be h^ -^ (X)i T^i)lJ'^ instead of h^ -* {Y.i ^i)^^- 

Fig.Elshows the correlation region for Oij ^ fij in the small O33 case. The 
curves correspond to the ratio Br{hf -^ Cl2i^i)f''~^)/B^i^i ~^ iJ2i^i)'^~^)^ 
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tan 023 

Figure 7: Br {hi -^ (^^ Vi)^^)/ Br{hl -^ (^^ i^i)T^){rn^/mrY as a function 
of tan 623 in the correlation region Oij ^ fij of the small O33 case. The solid 
curve corresponds to the best fit point value of tan 26^12, while the upper and 
lower curves corresponds to its 3a limits. See text. 

normalized by {m^/niT-y, as a function of the atmospheric mixing angle, as 
expected from Eq. (J53|l for the best fit point value of tan 29i2 (solid line) and 
its 3cr limits (dashed lines). The parameters are fixed as in Fig. IHland the 
spread of the points can be understood from the uncertainty in the solar mix- 
ing angle. In this region the charged Higss decay rate T{h~^ ^ {^- i'i)fi'^) can 
be larger than decay rate T{h'^ — ^ {J2i ^i)^~^) up to a factor of {171^-/171^)'^ = 
280. 

From Fig.iniit can be seen that the large O33 region is divided in three sub- 
regions, region I where 10~^ < Br{hf -^ {J2i ^i)l^^)/ Br{hl -^ (^^ ^i)'^^) ^ 
1, region II for which 1 < Br{ht ^ {T.i^i)^J^^)/ Br{ht -^ (Ei ^^^)^+) ^ 10^ 
and region III characterised by 10^ < Br{hf -^ {^i^i)^~^)/Br{hf —>■ 
(Ej ^i )''"''') ~ 10^- Measurementsof the ratio i?r(/;,]^ -^ {J2i^i)f''^)/Br{hf —^ 
(Ej ^i)''"^) ^^6 sufficient to decide whether region I or III are realized. In re- 
gion III there is an ambiguity that cannot be removed by measurements of 
the ratio Br{hf -^ (^i^^dt^^) / Brih^ -^ (Ej^j )''"''')• However, in the small 
mixing limit the ambiguity can be removed. Recalling that in the small O33 
region /12 = or /13 = one should expect, if this region is realized. 



5rf = Brl + Brl or Brl = Br': + Br^ 



(54) 



Any deviation from these relations would exclude this region and in addi- 
tion with a measurement of the type 1 < Briji^ —> C^ii^i)fJ'~^)/Br{h'^ — >• 
(Ej ^i)''''^) ~ 10^ "^ill indicate that region II is realized. 
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Figure 8: Br(h\ -^ {J2i^i)f^^)/B''"i^i ~^ Cl2i^i)'''^) ^s a function of O33 
and several pairs of Mi and M2 with fi = 100 GeV. From left to right the 
curves have sin2(y9 = 0.04, 0.15, 0.17 and 0.99. The dotted line is the same 
that the curve in the left part of Fig. IHl See text. 

The curves in Fig.lHlare basically independent of the value of fi. However, 
along each curve, smaller values of O33 are excluded as /i decreases. On the 
other hand they depend on the specific value of Mi and M2. In fact, as 
the mixing angle sin2(y9 increases the curves are shifted to the right. This is 
illustrated in Fig. |H1 for the small 0-^3 case with /12 = 0. All the remaining 
parameters are chosen as in Fig. IHl In particular the dotted line is the same 
as the one in Fig. IHl 

In summary for the one parameter-solutions we have 



m,. 



771, 






(55) 



We have checked that this result holds for all the parameter space of the 
NMZM. In particular for O33 sufficiently small the charged Higgs decay rate 
T{hf —* {J2i^i)f''^) '^^^ be dominant. 

7 Conclusions 

We have considered the version of the Zee model where both Higgs doublets 
couple to leptons. Instead of working with all the parameters we have fo- 
cused on a model with minimal number of couplings consistent with neutrino 
physics data. We have shown that in the small mixing limit {ip <^ 1) certain 
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ratios of branching ratios can be used to obtain information about the pa- 
rameters of the modeL Besides the charged Higgs leptonic decays we have 
also considered the decay /i^ -^ h^h^. We have found that this decay, if 
kinematically allowed, can be used to determine the value of the n parame- 
ter. Moreover, measurements of r(/i^ -^ h^h^) allow to decide whether the 
small mixing limit is realized or not. 

Assuming that there are no large hierarchies among the couplings On {i = 
1,2,3) and Oij, and using neutrino physics constraints we have shown that 
in this scheme only three parameters are independent. We have found that 
there are four regions, in this three-dimensional parameter space, determined 
by only O33. We have shown that two of these four regions are governed by 
large values of O33 (10"'^ — 10~^) while the other two regions are governed 
by small values of O33 (lO"^ - 10"^). 

We have analysed charged Higgs leptonic decays in the large as well as 
in the small O33 regimes and we have found: (i) in the large O33 case, there 
is a region in which the decays hf — >■ i^j/i"'" and hf -^ ViT^ are governed by 
the correponding Yukawas as in the 2HDM of type-I and type-II and another 
region where the decay h^ -^ Vi^Ji^ is enhanced and moreover can be larger 
than the decay to h'l —>■ ViT^ . (ii) In the small O33 case the decay h^ —>■ z/j/i+ 
is always enhanced and is larger than the decay hf —>■ ViT^ . Therefore we 
suggest that in order to test the model the decays of the charged Higgs to 
Viix^ should be searched along with the decays to ViT^ . In fact, measurements 
of the ratio of branching ratios Br{hf —>■ {^ii^i)f^'^)/Br{h^ -^ (Xli^*)''"^) 
could give information about what region of this parameter space is realized. 

At future colliders the decay channel z/r^ is very important for the dis- 
covery of charged Higgs bosons [211 ESI- For the LHC and SUSY like 2HDM, 
it has been claimed that the existence of a relatively heavy charged Higgs 
bosons, of mass up to 1 TeV, can be probed using the signal hf -^ z/r"*" 
|2l] . At future linear colliders a single produced charged Higgs should be 
associated with the tau and the neutrino coming from the virtual charged 
Higgs decay |2SI- According to our results, and illustrated by Fig. IHl the 
charged Higgs could emerge from a signal with z/j/i+ instead of I'iT'^. More- 
over, for a light charged scalar (Mi < rrit) the ratio of branching ratios, 
Br{hf -^ {J2i^i)f''~^)/^^i^i ~^ (Si'^i)'^^) should be measurable. 
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A The Three-parameter solution 

From the set of Eqs. p9|l we choose to express fu, /13, /23 and O13 in terms 
of O33, O22, and O12 
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